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'We present a study of various aspects of proximity effects in F/ S (Ferromagnet/Superconductor) 
bilayers, wliere F has a spiral magnetic texture such as that found in Holmium, Erbium and other 
materials, and S is a conventional s-wave superconductor. "We numerically solve the Bogoliubov-de 
Gennes (BdG) equations self-consistently and use the solutions to compute physical quantities rele- 
vant to the proximity effects in these bilayers. We obtain the relation between the superconducting 
transition temperature Tc and the thicknesses dF of the magnetic layer by solving the linearized 
BdG equations. We find that the Ta^dp) curves include multiple oscillations. Moreover, the system 
may be reentrant not only with dp, as is the case when the magnet is uniform, but also with temper- 
ature T: the superconductivity disappears in certain ranges of dp or T. The T reentrance reported 
here occurs when ciF is larger than the spatial period of the conical exchange field. We compute 
the condensation free energies and entropies from the full BdG equations and find the results are 
in agreement with Tc values obtained by linearization. The inhomogeneous nature of the magnet 
makes it possible for all odd triplet pairing components to be induced. We have investigated their 
properties and found that, as compared to the singlet amplitude, both the m = and m = ±1 
triplet components exhibit long range penetration. For nanoscale bilayers, the proximity lengths for 
both layers are also obtained. These lengths oscillate with dp and they are found to be long range 
on both sides. These results are shown to be consistent with recent experiments. We also calculate 
the reverse proximity effect described by the three dimensional local magnetization, and the local 
DOS, which reveals important energy resolved signatures associated with the proximity effects. 

PACS numbers: 74.45.-l-c, 74.62.-c, 74.25. Bt 



I. INTRODUCTION 

The emerging field of spintronics has stimulated in- 
terest in fabricating solid state devices that make use 
of the intrinsic spins as a degree of freedom. ^ Strides 
have been made recently towards exploiting the spin vari- 
able in hybrid ferromagnet/superconductor (F/S) sys- 
tems. Such systems have shovifn promise for a number of 
practical applications, including nonvolatile information 
storage. The simplest of such potential devices usually 
involve layered F/S heterostructures. Owing to these po- 
tentially important spintronic applications, research on 
the fundamental physics of these systems has received 
great attention in the last decade. ^'^ The most impor- 
tant basic physics elucidated by these studies is prob- 
ably that of the superconducting proximity effects in 
F/S nanostructures,'^ which describes the interplay be- 
tween ferromagnetic and superconducting order parame- 
ters. Though these two order parameters are rarely found 
to coexist naturally in bulk materials, such coexistence 
can be and has been achieved near the interfaces of ar- 
tificially created F/S composites. Thus, the subject has 
become important not only for its technological appli- 
cations but also because of the underlying fundamental 
physics. 

In elementary treatments, ferromagnetism is often 
deemed strictly incompatible with s-wave superconduc- 
tivity due to their mutually exclusive order parameters. 
In ferromagnets, the exchange field tends to cause the 
electronic spins to align in the same direction, while 



in singlet s-wavc superconductors, the Cooper pairs 
are composed of both spin-up and spin-down electrons. 
These two order parameters seem to naturally oppose 
each other. In fact, in F/S heterostructures this com- 
petition leads to a strong modification of the behavior 
of the superconducting Cooper pair amplitudes. When 
a Cooper pair in S encounters an F/ S interface and en- 
ters the F region, the momenta of spin-up and spin-down 
electrons that make up the Cooper pair are changed, be- 
cause of the exchange field in the F region. This leads to 
a nonzero center of mass momentum of the Cooper pair^'^ 
and an overall damped oscillating Cooper pair amplitude 
in the F side. It is because of these two competing order 
parameters that the oscillations decay over a relatively 
short length scale which decreases as the exchange field 
increases. 

These oscillations of the superconducting wavcfunc- 
tions are one of the most salient features governing prox- 
imity effects in F/S systems and form the basis for 
switching applications that require the manipulation of 
the superconducting transition temperature Tc through 
variation of the experimental parameters. Due to the 
oscillatory nature of the Cooper pair amplitudes, the 
dependence of Tc on the thickness of the ferromagnetic 
layer, dp, in F/S bilayers is oscillatory too. Furthermore, 
the interference between the transmitted pairing wave 
function through the F/S interface and the reflected one 
from the boundary can become fully destructive: the su- 
perconductivity disappears for a certain dp range. This 
superconducting reentrant behavior with dp has been 
found experimentally in Nb/Cui-^^Nia; and Fe/V/Fe 
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trilayers^® and it is well understood theoretically.'^'^ 

Another important fact about F/ S proximity effects 
is the generation of induced triplet pairing correlations. 
These can be generated by the presence of spin active 
interfaces, or (and this is the case we will focus on 
here) in systems with clean interfaces and inhomogeneous 
F structure^®"^^. The simplest such cases are F1/S/F2 
or F1/F2/S layers in which the magnetizations of the 
two F layers are misaligned. For s-wave superconduc- 
tors, where the orbital part of the pair wave function is 
even, the Pauli principle requires the spin part to be odd 
and this would appear to forbid the existence of triplet 
correlations. However, triplet correlations that are odd^^ 
in frequency (or equivalently in time^^) can be induced 
in F/S systems, with S being s-wave pairing, without 
violating the Pauli exclusion principle. 

The importance of odd triplet correlations lies in their 
long range nature in the magnet, i.e., their proximity 
lengths can be in principle comparable to those found 
in the usual superconducting proximity effects involving 
nonmagnetic metals. Since the exchange fields tend to 
align the electronic spins of the Cooper pair electrons, 
the proximity length for singlet pairing is very short (and 
dependent on the magnitude of exchange field) . However, 
the triplet pairing correlations can involve electron pairs 
with both spins aligned along the local magnetization 
direction, and thus be much less sensitive to the mech- 
anism of exchange fields, penetrating much deeper in F 
than their singlet counterparts. The possible appearance 
of both m = and m = ±1 (m denoting the usual spin 
quantum number) components of induced triplet corre- 
lations is controlled by the symmetry of the system and 
by conservation laws^^. In multilayer F/S systems, when 
the F layers are magnetically homogeneous (all exchange 
fields are along the same direction, the quantization axis) 
only the total spin projection corresponding to the m = 
component can be induced. On the contrary, all three 
components (m = and m — ±1) can arise if the direc- 
tion of exchange fields differs in the ferromagnets, e.g. the 
exchange fields of Fi and F2 are not aligned in^^ Fi /S/ F2 
or F1/F2/S types of trilayers^*. These long-range charac- 
teristics of triplet correlations have been experimentally 
detected in ferromagnetic multilayers by taking advan- 
tage of their magnetic inhomogencity.^^"^^ 

Besides ferromagnet misalignment, another possibil- 
ity to generate long-range triplet correlations is to use 
a ferromagnet with an intrinsic inhomogeneous mag- 
netic texture.^® Such structures are inherent to cither 
known elements or chemical compounds. Examples of 
this kind of ferromagnets include most prominently Ho,^^ 
which has a spiral magnetic structure at low temper- 
atures. A similar spiral magnetic structure is found 
in metallic Erbium^°, MnSi thin films^\ and Fe(Sc,Tc) 
compounds'^^. Indeed it has been experimentally con- 
firmed that long-range triplet correlations are induced in 
Nb/Ho/Co multilayers'^'^ with the periodicity of Ho play- 
ing an important role in triplet supercurrcnts. Super- 
conducting phase-periodic conductance oscillations have 



also been observed in Al/Ho bilayers^^ where the thick- 
ness of Ho is much larger than the penetration length of 
singlet amplitudes. This finding can be explained in the 
framework of the triplet proximity effects. Theoretically, 
the spin-polarized Josephson current in S/Ho/S junctions 
has been studied^^ via quasi-classical Green function 
techniques. The triplet supcrcurrent in Ho/Co/Ho tri- 
laycrs was also investigated in the diffusive-^^ and clean'^^ 
regimes. The long-range effects can however be limited 
by interface quality and impurities.'^'' These earlier works 
show that ferromagnets with an intrinsic spiral magnetic 
structure are of particular interest when studying super- 
conducting proximity effects in F/S nanostructures. 

It was also recently predicted'^^ that superconductivity 
in conical- fcrromagnet/superconductor bilayers can be 
reentrant with temperature, in addition to the standard 
reentrance with dp mentioned above. It was shown via 
numerical solution of the Bogoliubov-de Gennes (BdG) 
equations that in certain cases superconductivity can ex- 
ist in a range Td < T < Tc2, where Td is nonzero. This 
reentrance with temperature is quite different from that 
found in ternary rare earth compounds such as ErRh4B4 
and HoMoeSs^^"'*'^ where the disappearance of supercon- 
ductivity below Tci results from the onset of long-range 
ferromagnetism. In the bilayers considered in Ref. 38 
the high T phase and the low T phase are the same. The 
physical reasons that account for the reentrance there are 
attributed to the proximity effects associated with the 
interference of Cooper pair amplitudes and the genera- 
tion of triplet pairing correlations, resulting in nontrivial 
competition between the entropies and condensation en- 
ergies. 

In this paper, we present results for various properties 
of the proximity effects in F/S bilayers, where the F layer 
has a helical magnetic structure. We numerically find the 
self-consistent solutions to the BdG equations'*^ and use 
them to compute important physical quantities. By lin- 
earizing the BdG equations we calculate the critical tem- 
perature as a function of magnet thickness, exchange field 
strength and periodicity, and other parameters. We then 
discuss the effects of varying the superconductor thick- 
ness to coherence length ratio. We show that depending 
on the width of the superconductor, and for a broad range 
of magnetic strengths, reentrant behavior as a function 
of magnet thickness can arise. We find that under certain 
conditions, the superconductivity can also be reentrant 
with temperature, and for larger dp values than previ- 
ously reported.'^® To clarify these reentrant phenomena, 
we investigate the thermodynamic functions associated 
with the various ways reentrance can arise. We find that 
all components of the odd triplet correlations can be in- 
duced and wc discuss their long range nature. We then 
characterize the important triplet long range behavior 
by introducing the corresponding proximity lengths. We 
find that these lengths oscillate as a function oi dp, and 
depend on details of the magnetic texture. Reverse prox- 
imity effects are also studied to determine the magnetic 
influence on the superconductor: we calculate the local 
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FIG. 1. (Color online) Diagram of the conical ferromagnet- 
superconductor bilayer studied. The spiral magnetic struc- 
ture is described by an exchange field h (see Eq.(l)). The 
system is infinite in the x — z plane and finite in y. The 
relevant widths are labeled. 



magnetization vector revealing greater penetration into 
S for weaker exchange fields. Lastly, the spectroscopic 
information is presented by means of the local density of 
states, to demonstrate consistency with the results. 



II. METHODS 

The procedures we employ to sclf-consistcntly solve the 
BdG equations and to extract the relevant quantities are 
very similar to those already described in the literature 
(see Refs. 17 and 38 and references therein). It is unnec- 
essary to repeat details here. We consider F/S bilayers 
that consist of one ferromagnetic layer with spiral ex- 
change fields and a superconducting layer with s-wave 
pairing. The geometry is depicted in Fig. 1. Our sys- 
tems are infinite in the x — z plane and finite along y-axis 
Their total thickness is denoted by d: the F layer has 
width dp and the S layer has width ds = d — dp- The 
left end of the bilayers is the y = plane. We assume 
that the interface lies in the x — z plane and the exchange 



field h, which is present only in F, has a component that 
rotates in this plane plus a constant component in the y 
direction, perpendicular to the interface: 



h. = ho I cos ay -I- sin a 



sm I — X + cos — I z 

ay V a , . 

(1) 

where the helical magnetic structure has a turning an- 
gle /3, and opening angle a. We will take a, the lattice 
constant, as our unit of length and vary the strength Hq. 
The spatial period of the helix is A = 27ra//3. 

We first write down the effective Hamiltonian of our 
systems. 
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where p and p' are spin indices, V'p(r) is the field oper- 
ator, and (T are the Pauli matrices. A(r) in the above 
equation is the pair potential. To apply the BdG formal- 
ism to spatially inhomogeneous systems, we first invoke 
the generalized Bogoliubov transformation. 



(3) 



where u„p(r) and f„p(r) can be interpreted as quasi- 
particle and quasi-hole wavefunctions and the creation 
operator and annihilation operator 7„ obey the usual 
fermionic anti-commutation relations. To recast the ef- 
fective Hamiltonian into a diagonalized form, Heff = 
J2n ^nlnln; ouc Can usc Commutation relations between 
Heff and field operators. By doing so, and making use 
of the quasi one dimensional nature of the problem one 
arrives at the BdG equations. 
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where Ho = — ^i—^ ■§jj7 + e± — Ep , is the usual single par- 
ticle Hamiltonian for the quasi one-dimensional problem, 
with e_L denoting the kinetic energy associated with the 
transverse direction. By using Eq. (3), the self consis- 
tency relation for the pair potential can be rewritten in 



the form, 

A(y) = Kt(yX4(y) ~ ^ni{y)vntiy)] tanh(|^), 

n 

(5) 

where g{y) is the usual BCS superconducting coupling 
constant in the S region, and zero in the F material. The 
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prime sign indicates summing over all eigenstates with 
eigenenergies less than or equal to a cutoff "Debye" fre- 
quency uji) . The solutions to the BdG equations must be 
determined self-consistently. This self-consistency condi- 
tion is extremely important in studying proximity effects. 

The singlet pair amplitude, i.e., the amplitude for find- 
ing a Cooper pair, F{y), is given by F{y) = A{y)/g{y). 
One can determine the superconducting transition tem- 
peratures Tc by looking for the temperatures at which the 
pair amplitudes becomes vanishingly small. However, it 
is much more efficient to find Tc by linearizing the self 
consistency relation, Eq. (5) and using a perturbation 
expansion. This technique has been discussed in other 
papers-^®'"*^ and the details are not repeated here. 

Once a full set of self-consistent solutions is obtained, 
all the additional quantities of interest can be computed. 
For example, the triplet correlations, can in general be 
defined^^'^^ as: 

/o(r,i) = i[(^t(r,OV'i(r,0)) + (^;(r, t)Vt(r, 0))], (6a) 

A(r,t) ^ i[(^^(r,<)^t(r,0)) - (^;(r, i)V'i(r, 0))], (6b) 

where (...) represents ensemble averages. As discussed 
in Sec. I, both /o and /i have to vanish at t = to com- 
ply with the Pauli principle. However, in general both of 
them can be induced when t 7^ and the magnetic struc- 
ture is inhomogeneous. By using Eq. 3 and considering 
the time evolution, one can rewrite the odd triplet cor- 
relations in terms of quasi-particle and quasi-hole wave- 
functions in the form: 

n 

(7a) 

h{y,t) = [""t(y)<t(y) -""4.(yXi(y)] Cn{t), 

n 

(7b) 

where Cn{t) = cos(e„i) — i sin(e„t) tanh(e„/(2r)). 

Given the quasi-particle amplitudes and eigenvalues, 
one is also able to evaluate the thermodynamic quantities 
from the free energy F{T). For an inhomogeneous system 
it is most convenient to use the expression:^^ 

f(T).-2TX;i.[2co.h(i^)]+(^)_, (8) 

where here (. . .)s denotes spatial average. The conden- 
sation free energy AF is the difference between the free 
energies of the superconducting state, Fs, and the nor- 
mal states Fn, i.e. AF = Fs — F^. Fn can be calculated 
by assuming that the pair potential is absent throughout 
the system. 

Another important physical quantity, which can be 
determined experimentally by tunneling spectroscopy, is 
the local density of states (LDOS). This quantity often 



reveals important information about the superconduct- 
ing features of the sample studied. In our quasi-one- 
dimensional model, the LDOS N{y, e) depends spatially 
only on y. N{y, e) consists of both spin-up and spin- 
down LDOS contributions, that is, N{y,e) = N-f{y,e) + 
Ni{y,e). 

Npiy, e) = ^[|u„p(2/)p5(e - e„) + \vnp{y)\^Sie + e„)], 

n 

(9) 

where p =t, l- 

Just as the superconducting order parameter is 
changed by the presence of ferromagnets, also, near 
the interface, the ferromagnetism can be modified by 
the presence of the superconductor, a phenomenon 
known as the reverse proximity effect. It is best described 
by considering the local magnetization m. The local 
magnetization vector is defined as m = — Mb(X](t ^'^cr^'), 
where ^ = {tpi^,ipi)'^ , and, again, it depends on the co- 
ordinate y only. By using Eq. 3 we have, 

rrixiy) = " ^-B X!"^ (Ktiyhniiy) + <^(y)"nt(2/))/n 

+ iv,it{y)vni{y) + vni{y)vlf{y)){i - /„)|, 

(10a) 

my{y) ==i^B ^1 {ul^{y)uni{y) - Kiivhuf {y))fn 

+ {Vnt{y)Vnl{y) - Vnl{y)v^^{y)){l - /„)|, 

(10b) 

mziy) = - A^B X! 1 (l"»t(y)P - Wnliy)\'^)fn 
n ^ 

+ (Kt(2/)l'-k4(y)l')(i-/n)}, (loc) 

where /„ is the Fermi function of e„ and /^s is the Bohr 
magneton. 



III. RESULTS 

In the results shown here, capital letters will always 
denote the dimensionless lengths denoted by the corre- 
sponding small letter. For example, the dimensionless 
thickness of the fcrromagnct is written as Dp = dp /a 
and that of superconductors is Ds = ds/a, where a is 
the lattice constant in Eq. (1). For the helical magnetic 
structure we take angular values (see Eq. (1)) a = 47r/9 
and /3 = 7r/6 which are^^'^^ appropriate to Ho, in which 
case Dp = 12 contains one full period of the spiral ex- 
change field. We will denote this dimensionless spatial 
period by A in the following subsections. For materi- 
als other than Ho many of the results can be read off by 
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FIG. 2. (Color online) Calculated transition temperatures T^, 
normalized to r°, vs. Dp for several values of the dimension- 
less exchange field I (see text). In this figure Ds is fixed for 
all values of I to be f .5Ho. The lines connecting data points 
are guides to the eye. 
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FIG. 3. (Color online) Transition temperatures Tc vs. Dp at 
/ = O.I and several Ds- The lines are guides to the eye. 



rescaling A to the appropriate value. Throughout this pa- 
per, the dimensionless superconducting coherence length 
is fixed to be Sq = 100. In the same spirit, the dimension- 
less exchange field, /, is measured in terms of the Fermi 
energy: I = ho/ Ep ■ We choose the Fermi wavevector in 
S to equal 1/a. Wc take the "Debyc" cutoff value to be 
lod = Q.QAEp- As usual, this value is irrelevant except 
for setting the overall transition temperature. Temper- 
atures are given in dimensionless form in terms of T^, 
the transition temperature of bulk 5" material. When 
discussing the triplet amplitudes, which are time depen- 
dent, we use the dimensionless time t = tupi. Vertical 
dashed lines shown in figures, when present, denote the 
F/S interface. 



A. Transition Temperatures 

To investigate the details of the predicted oscillatory 
nature of the dp dependence of as discussed in Sec. I, 
we calculated as a function of Dp for several / and 
Ds values. These results are shown in Fig. 2 and Fig. 3. 
The Dp range in both figures includes three complete 
periods of the spiral magnetic order. This is reflected in 
the results shown: indeed, the presence of multiple oscil- 
lations in the included range of I? is the most promi- 
nent feature in Figs. 2 and 3. The oscillations in Tc arise 
(as we discuss below) from a combination of the period- 
icity of the spiral magnetic structure and the usual 
oscillations which arise, even when the magnet is uni- 
form, from the difference'^'^ in the wavevectors of the up 
and down spins. In Fig. 2, one can also see that with 
stronger exchange fields the oscillation amplitudes are 
larger. Despite this increase of the amplitudes with the 
exchange field (they are approximately proportional to 
/), the overall Tc decreases when the exchange field in- 
creases. This is consistent with expectations: a stronger 
exchange field destroys the superconductivity more effi- 
ciently. As mentioned in Sec. I, when the exchange field is 
strong enough, the systems can become normal in some 
range Dpi < Dp < Dp2. Indeed, reentrance with dp 
can be seen to occur in Fig. 2 near Dp — i a,t I ~ 0.2. 
Another feature seen in this figure is the decrease of the 
amplitude oscillations with increasing Dp. This arises 
simply because the singlet Cooper pair amplitudes in S 
near the F/ S boundary decay more strongly at larger dp 
and therefore the effect of the pair amplitude oscillations 
in F is weaker. 

In a F/S bilayer where the ferromagnet is homoge- 
neous, the periodicity of the Tc oscillations is governed 
by the exchange field, or equivalently, by the magnetic co- 
herence length^ Ep ~ l/I. Here, where a bilayer with a 
conical inhomogeneous ferromagnet is considered, the in- 
trinsic spiral magnetic order with spatial period A plays 
an equally important and competing role in the Tc os- 
cillations. In other words, both the strength and the 
periodicity of exchange fields influence the overall decay 
and the oscillatory nature of the superconducting transi- 
tion temperatures. The existence of two different spatial 
periodicities leads to the obvious consequence that the 
Tc{Dp) curves are not describable in terms of one single 
period. However, when / is not very strong (/ < 0.1), 
the minima of Tc are near the locations where Dp = A/2, 
3A/2, and 5A/2 and similarly, the Tc maxima occur near 
Dp = A, 2 A, and 3 A. This indicates that the magnetic 
periodicity is dominant. Roughly speaking, the maxima 
and the minima are correlated with the strongest and 
weakest spatial average of the exchange field components 
in F. As / increases and Sj? decreases deviations become 
obvious. Figure 2 shows that the distances between two 
successive maxima decrease when the exchange fields in- 
crease. 

The existence of the multiple oscillations discussed 
above has been confirmed experimentally. In Ref. 53, 
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FIG. 4. (Color online) Calculated transition temperatures Tc 
vs. Dp for Ds ~ 148 and / — 0.2. The main plot shows ((red) 
symbols) the overall Tc behavior from Df = Q to Df = 1.5A. 
Reentrance with Df near Df = 4 is seen. In this case there 
is also reentrance with temperature in the region indicated 
by (green) symbols near Df~W- The inset is a blow up 
of this region: superconductivity exists only in the region 
Tci < T < Tc2, where Tc2 is depicted by the upper (green) 
squares and Tci by the (blue) circles. 



Tc in Nb/Ho bilaycrs was measured as a funetion of dp- 
The results exhibit an overall decay with Ho thickness, 
on which there are superimposed oscillations which arc 
correlated with, but not simply described by, the spatial 
wavelength A of the Ho structure. Comparison with the 
theory discussed here was made, using / as an adjustable 
parameter. Values near / = 0.1 were found to provide the 
best fit. The other parameters were extracted from other 
known properties of Ho and Nb or (e.g. ds) from the ex- 
perimental sample geometry. The results of the compari- 
son were extremely satisfactory, showing clear agreement 
in all the features of the rather intricate Tc{dF) experi- 
mental curves. It was found also that one of the samples 
was close to being reentrant with c?i^ at a value very close 
to that predicted by theory. 

In Fig. 3, we present Tc results for several values of 
Ds ranging from 1.2So to 1.5So, with a fixed exchange 
field / = 0.1. One can see that the distance between 
successive maxima is an extremely weak function oi Ds- 
This agrees with our previous discussion: the oscillatory 
nature in is chiefly dependent on the exchange fields 
and magnet structure. Since superconductivity is more 
robust for larger Ds, the ferromagnet lowers the over- 
all Tc for thinner superconductors as evidenced in Fig. 3. 
Fig. 3 also demonstrates that not only a strong / can lead 
to Dp reentrances, but also a thinner Ds- Interestingly, 
at the smallest value of Ds considered, there are two Dp 
reentrance regions, one near Dp = 5 and the other near 
Dp = 27. As discussed above, these Dp reentrances in 
both Fig. 1 and 2 are mainly due to the interference 
between the transmitted and reflected Cooper pair con- 



densates that are oscillatory in the F region. 

In previous work^®, we reported one speciflc case where 
superconductivity in Ho/S bilayers exhibits not only the 
usual reentrance with dp but also, at some flxed values of 
ds, ho, and dp, reentrance with T, that is, superconduc- 
tivity exists only in a temperature range Td < T < Tc2, 
where Td is flnite. In the example reported in Ref. 38 
temperature reentrance occurred near the flrst minimum 
of the Tc{Dp) curve. We have investigated here whether 
this kind of reentrance can occur near some of the other 
minimum of Tc{Dp)- These locations appear favorable 
for such an occurrence since superconductivity is rela- 
tively weak near these minima. Also reentrance with Dp 
is after all an extreme case of a minimum Tc{Dp)- We 
have found that other T-reentrant examples can indeed 
be found, although by no means universally. Here we 
report an example of reentrance occurring near the sec- 
ond minimum of Tc{Dp)- In Fig. 4, the main plot shows 
Tc{Dp) for the parameter values spcciflcd in the caption. 
The flrst minimum of Tc{Dp) drops to zero and is an 
example of Dp reentrance. In the region near the second 
minimum ((green) symbols) reentrance with T occurs. 
The region of interest is enlarged in the inset. There 
the upper (green) symbols represent Tc2 and the small 
dome of lower (blue) circles represent Td- In the dome 
region, but not outside it, the superconductivity is reen- 
trant in T - When one lowers the temperature from the 
normal region, the F/ S bilayers become superconducting 
at Tc2- With further cooling, the bilayers return to nor- 
mal state. Reentrance in this case occurs at the second 
minimum rather than the first because there is no upper 
transition associated with the first minimum: the system 
is normal. Near the second minimum the oscillatory ef- 
fects are not as strong, and as as a result, the system 
becomes reentrant in T. This can be viewed as a "com- 
promise" : near the second minimum, as opposed to the 
first, superconductivity is not completely destroyed but 
it becomes "fragile" and can disappear upon lowering T - 
The physics involved from the thermodynamic point of 
view will be discussed in the following subsection. 



B. Thermodynamics of Reentrance phenomena 

To understand the reentrance phenomena in T it is 
most useful to examine the thermodynamics of the two 
transitions, and in the region between them. From the 
condensation free energy A_F, which can be evaluated 
as explained in connection with Eq. (8), other quantities 
such as the condensation energy and entropy are easily 
obtained. For reentrance with Dp, it is sufficient to look 
at the free energy at constant low T. 

Considering the reentrance with T , it is illuminating 
to consider first the T dependence of the singlet pair 
amplitude F{Y) well inside the S material. Thus, we 
focus on F{Y) one coherence length from the interface: 
Y = Dp + 'E.Q- This quantity, normalized to its value in 
bulk S material, is plotted in Fig. 5 as a function of T 
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FIG. 5. (Color online) The singlet pair amplitude, normalized 
to its value for bulk S material, at a location one coherence 
length inside S from the Ho/S interface, plotted as a function 
of T. The (red) squares are for Dp = 16 and the (green) 
circles are for Dp = 17. All other parameters are as in Fig. 4. 



for two contrasting values oi Dp, one at Df = 16 where 
rccntrancc occurs (sec Fig. 4) and at a very nearby value, 
Dp = 17, which lies just outside the reentrance region 
and exhibits typical behavior. Fig. 5 ((green) circles) 
demonstrates that in the latter case the amplitude be- 
haves qualitatively as the order parameter does in a con- 
ventional BCS superconductor: it decreases very slowly 
near T ~ and eventually drops to zero very quickly 
but continuously near Tc, indicating the occurrence of a 
second order phase transition. This transition occurs at 
Tc/T^ ~ 0.32 in agreement with Fig. 4. However, the 
behavior of the pair amplitude in the reentrant region 
((red) squares in Fig. 5) is quite different. There are 
two transition temperatures: below a very low but finite 
temperature, T^i/T^ = 0.02, the singlet pair amplitude 
vanishes and the system is in its normal state. F{Y) 
then begins to rise continuously, has a maximum at a 
temperature T^, (where Tm/T^ ~ 0.1) and eventually 
drops to zero, again continuously, at an upper transition 
Tc2/T° « 0.22. In the region Tci<T < Tc2 the system is 
in the superconducting state. Both transitions are of sec- 
ond order. The values of Td and Tc2 from the vanishing 
of the amplitude, seen in Fig. 5, agree with those cal- 
culated directly from linearization of the self consistent 
equation plotted in Fig. 4. 

We now turn to the condensation free energy, AF, and 
entropy, AS*, for the same T reentrant case. AF is shown 
in the top panel of Fig. 6 as calculated from Eq. (8) 
and normalized by 2£'o, where Eo is the condensation 
energy of bulk S material at T = 0. The lower panel 
shows the normalized condensation entropy, defined as 
the AS = -dAF/diT/T^). The meaning of the symbols 
in this figure is the same as in the previous one. When 
the system is near (but outside) the reentrant region the 
behavior of both quantities plotted is qualitatively the 
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FIG. 6. (Color online) The normalized condensation free en- 
ergies, A_F = Fs — Fn, vs. T/T^ are shown in the top panel 
for the same cases presented in Fig. 5. The (red) squares and 
right scale are for Dp — 16. The (green) circles and left scale 
are for Dp ~ 17. The bottom panel shows the normalized 
(see text) entropy differences, AS = Ss — Sn vs. T/Tj^, on 
the same vertical scale. The meaning of the symbols is the 
same as in the top panel. 



same as that found in textbooks for bulk BCS supercon- 
ductors. Quantitatively, the magnitude of AF for our 
systems are much smaller than that for bulk S where we 
would have AF = — 0.5atT = 0in our units. The value 
of Tc in the non-reentrant case can also be identified from 
where the free energies of the normal and superconduct- 
ing states are the same {AF{T) = 0), and it agrees with 
both Fig. 4 and 5. Moreover, the vanishing of the en- 
tropy difference at a finite Tc confirms the occurrence of 
a second order phase transition. The value of this tran- 
sition temperature is consistent with all above results. 

The story for the reentrant case is quite different. 
There, although the values of AF are much smaller com- 
pared to those in the standard case, one can still find 
that the minimum of AF occurs at approximately the 
same value T,„ where the singlet pair amplitudes have a 
maximum. Thus the superconductivity is most robust at 
T = Tm- The two transition temperatures Td and Tc2 
can also be determined from the top panel of Fig. 6 and 
match with those found in Figs. 4 and 5. In the two 
T ranges T < Td and T > Tc2; the normal state is the 
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FIG. 7. (Color online) Reentrance with Df- Top panel: nor- 
malized singlet pair amplitude, computed at a location one 
coherence length inside S from the interface, as a function 
of _Df. Bottom panel: normalized condensation free energy, 
AF = Fs ~ Fn, vs. Df at T = 0. 



only self-consistent solution to the basic equations, as is 
evident from Fig. 5. The vanishing AF when T < Td 
means that the electrons do not then condensate into 
Cooper pairs. This is exactly what happens for pure su- 
perconductors when T > Tc- 

There are some remarkable facts about the behavior 
of AS" in the reentrant case. First, the vanishing of AS' 
(along with that of AF) in Fig. 6 indicates that the sys- 
tem undergoes second order phase transitions at both Td 
and Tc2- Also, AS" is positive for Td < T < Tm where 
is again the value of T at which the singlet pair am- 
plitude reaches its maximum and AF its minimum. That 
the entropy of the superconducting state is higher than 
that of the normal state indicates that the normal state 
at Tci < T < Tm is more ordered than the superconduct- 
ing one. This truly unusual fact, which is the root cause 
of the reentrance, is due to the oscillating nature of both 
the Cooper pair condensates and of the exchange field, 
which leads to an uncommonly complicated structure for 
the pair amplitude. Above Tm, the superconducting state 
becomes more ordered than the normal state: AS is neg- 
ative. From Fig. 5 and 6, we see that the singlet pair 
amplitudes, the condensation free energies, and entropy 
differences of reentrant case in the range Tm < T < Tc2 



have a similar trend to those of the non-reentrant case in 
the range < T < T^. We have found also examples of 
non-reentrant cases in which there is a finite temperature 
Tm at which AF has a minimum but where on further 
lowering T , AF remains negative all the way to T = 0. 

The situation in the more common Dp reentrance re- 
gion, where we find that the system does not become 
superconducting when it is heated from T = 0, is differ- 
ent from that of T reentrance. A case where Tc vanishes 
in the range Dpi < Dp < Dp2 ioi I = 0.2 was seen in 
Fig. 2. To further analyze this Dp reentrance, we again 
calculated the singlet pair amplitudes inside 5* at one 
coherence length from the interface, in the zero tempera- 
ture limit. The top panel of Fig. 7 shows the normalized 
F{Dp + 'E.q) as a function of Dp, for the same parameters 
as the I = 0.2 case in Fig. 2. The singlet amplitudes drop 
to zero in the same range as where T^ vanishes in Fig. 2: 
the normal state is the only self-consistent solution and 
the superconductivity is completely destroyed in this D p 
range. One can also see that the order parameter is con- 
tinuous but its derivative discontinuous at Dpi and Dp2. 
In the bottom panel, we plot the corresponding conden- 
sation free energies (at T = 0) as a function oi Dp. The 
Dp range and the temperature are the same as the top 
panel. The condensation free energies vanish in the same 
Dp reentrance region. Unlike the derivatives of the sin- 
glet pair amplitudes, the derivatives of A_F at Dpi and 
Dp2 appear to be continuous. 

The physical origins of these two kinds of reentrance 
are not identical. As mentioned in Sec. I, the interference 
effects of oscillating Cooper pair wavefunctions are re- 
sponsible for the Dp reentrance, provided that / is strong 
and Dp is not too thick. Dp reentrance does not require 
a nonuniform magnet. The conical-ferromagnet struc- 
ture introduces an additional nonuniform magnetic order 
which may coexist with nonuniform superconductivity 
as predicted in Ref. 54. This additional non-uniformity, 
with its concomitant introduction of triplet correlations 
and of a new periodicity, can produce, as we have shown, 
reentrant behavior in T , as opposed to the simpler be- 
havior seen e.g. near the first minimum in the main plot 
of Fig. 4. Thermodynamically, the reentrance with T is 
due to the competition between entropy and energy,^® 
and driven by the high entropy of the disordered super- 
conducting state. When T < Tm, AS is positive and the 
roles of the normal and superconducting phases are ex- 
changed: the high entropy phase is the superconducting 
one. Further lowering T brings the system back to nor- 
mal state. One can compare the instance of T reentrance 
reported here with that reported in our previous work'^* 
where it occurs near the first minimum of Tc{Dp). In 
that work, Ds = 150 = 1.5So and / ~ 0.15. Here, not 
only is Ds thinner but also / is greater. The first mini- 
mum of Tc{Dp) in the main plot of Fig. 4 drops to zero 
and becomes & Dp reentrance region. Because stronger 
/ and thinner Dg are unfavorable to superconductivity, 
the system can not sustain T reentrance there. Thus, a 
delicate balance of geometrical and material parameters 
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FIG. 8. (Color online) Comparison between the spatial dependencies of the singlet pair amplitude and the induced triplet 
correlations at the two indicated values of I, at Ds = 150, Dp — 300 and T = 0. The S region is to the right of the 
dashed vertical line. Both singlet, F{Y) ((green) curves higher in the S region)) and triplet, fo{Y), fi{Y), pair amplitudes are 
normalized to the value of F{Y) in pure bulk 5* material. For this comparison, the normalized induced triplet pair amplitudes, 
which are evaluated at r = 9.6, are multiplied by a factor of 10. The real parts of fo{Y) and fi{Y) are shown ((red) curves 
strongly oscillating in the F region)). 



is required. 



C. Singlet to triplet conversion 

In this subsection, we will discuss the general proper- 
ties of the induced triplet pairing correlations in F/ S bi- 
layers with F being a conical fcrromagnet. As mentioned 
in Sec. I, in the presence of inhomogeneous exchange 
fields in the F layers both the to = and the m = ±1 
triplet pair amplitudes are allowed by conservation laws 
and the Pauli principle, but this says nothing about their 
size or shape, or indeed on whether they will exist at 
all. Thus detailed calculations are needed. The intrin- 
sically inhomogeneous magnetic textures discussed here 
provide unique opportunities to study the triplet prox- 
imity effects in F/S systems containing only a single F 
layer. Triplet correlations in the ballistic regime for both 
F1/S/F2 and F1/F2/S trilayers have been found in previ- 
ous work^^'^^'^* to be long ranged and the expectation'^'^ 
that they will also be in our case is fulfilled. We will here 



discuss and characterize this and other aspects (such as 
the effect of the strength of the exchange fields on the 
triplet pair amplitudes) of triplet pairing correlations in 
F/S bilaycrs where the magnets maintain a spiral ex- 
change field. Results presented in this subsection are all 
in the low T limit. 

To exhibit the long range nature of both types of triplet 
amplitudes, we show in Fig. 8 both the triplet and sin- 
glet pair amplitudes for a thick F layer as a function of 
position, as given by the dimensionless coordinate Y. In 
this (and the next figure. Fig. 9), we will focus on the 
real parts of the in general complex (see Eqs. (7)) /o 
and /i, since we have found that, for the cases shown, 
their imaginary parts are smaller by at least a factor of 
2 to 5 and their behavior is similar to that of the real 
parts. To properly compare singlet and triplet quanti- 
ties, both the singlet amplitude, F{Y) and the triplet 
amplitudes are normalized the same way: to the value 
of the singlet amplitude in bulk S material. For visi- 
bility, we have multiplied the triplet pair amplitudes by 
a factor of 10. In the left and right columns of Fig. 8, 
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FIG. 9. (Color online) Real parts of the induced triplet pair 
amplitudes, normalized as in the previous figure, for different 
characteristic times r. In these plots, Ds = 150, Dp = 300 
and / = 0.5. The top panel shows the real part of fo{Y) and 
the bottom one that of fi{Y). 



we show in this way the real parts of both /o and /i 
when / = 0.1 and / = 0.5, respectively. The ferromagnet 
has a large thickness: Dp = 3So = 20 s- The triplet 
correlations, which we recall vanish at equal times, are 
computed at a value of the dimensionless time t = 9.6. 
One sees right away that both the /q and /i components 
can be induced at the same time. This is always the case 
in our structures, as opposed to what occurs Ya.F1/F2/S 
and F1/S/F2 trilayers where the /o and /i components 
can be induced simultaneously only when the exchange 
fields in these F layers are non-coUinear. Secondly, the 
induced triplet correlations on the F side are long ranged 
compared to the singlet amplitudes. The singlet ampli- 
tudes decay with a short^ proximity length 27r!Bi? « 27r/J 
due to the pair-breaking effect of the exchange field. In 
contrast, the proximity length for the triplet amplitudes 
as seen in Fig. 8 is much longer: it is of the order of Sq, 
and docs not depend strongly on /. The triplet ampli- 
tudes spread over the F side with an oscillatory behavior. 
This difi^erence is more pronounced in the / = 0.5 case, 
where the decay length Sj? is much shorter than and 
the singlet amplitudes diminishes much faster than in the 
/ ~ 0.1 case. For both / = 0.1 and / = 0.5, one can also 
see that the singlet amplitudes begin to rise from the 



F/S interface and saturate in the S side about one su- 
perconducting coherence length from the interface. This 
agrees with our previous work^. Another interesting fea- 
ture seen in the / = 0.5 case is that the peak height of 
the /i component near the interface is not much higher 
than that of its other peaks, as happens with its fo coun- 
terpart. In other words, the subsequent peak heights in 
the F regions are comparable to that of the peak nearest 
to the interface. 

In delineating the role of triplet correlations in other 
experimentally relevant quantities, it is necessary to un- 
derstand their time dependence. Due to the self con- 
sistent nature of the proximity effects and the fact that 
the triplet condensate amplitudes are odd in time, their 
time dependence is in general nontrivial. We illustrate 
this in Fig. 9, where we show the spatial dependence of 
both the TO = and to = ±1 components of the triplet 
amplitude for several r. The parameters used here are 
the same as in the right panels (/ = 0.5) of Fig. 8. For 
easier comparison with Fig. 8, we have again multiplied 
the normalized triplet amplitudes by a factor of 10. Fig- 
ure 9 shows that at small times triplet correlations are 
generated only near the interface. (We have of course 
verified that they always vanish when r = 0). One can 
extract information about the proximity length from the 
growing increase of peak heights in the F regions. The 
peak heights grow faster when they are deeper inside the 
ferromagnet. Moreover, Fig. 9 clearly demonstrates that 
the triplet correlations penetrate into F regions as r in- 
creases, in the range studied. More remarkably, the peaks 
of the /i component that are not nearest to the interface 
grow very fast in time and have heights that are compara- 
ble to the one nearest to the interface, consistent with our 
remarks in our discussion of Fig. 8. In contradistinction 
with the oscillating behavior of the triplet amplitudes in 
the F regions, one can see that both fo and /i decay 
monotonically into the S side without any oscillations. 
However, the triplet correlations still spread over in the 
S regions at larger values of r just as they do in the F 
layer. 

In the above paragraphs, we have discussed the long 
range nature and other properties of the triplet ampli- 
tudes in our system when the conical ferromagnet is very 
thick. In the following paragraphs, we will consider the 
proximity effect of induced triplet pairing correlations for 
smaller scale conical- ferromagnets. To quantify the effect 
we introduce a set of proximity lengths Li^M defined as: 



L 



i.M 



J,jdY\MY,r)\ 
maxM |/i(y,r)| 



i = 0, 1 M^S,F. (11) 



Here the first index denotes the spin component, and 
the second index M denotes the region in which the 
given function is evaluated. If the decays were exponen- 
tial, these lengths would coincide with the characteristic 
length in the exponent. Obviously, in the present situa- 
tion the decays are more complicated but the Li M can 
easily be extracted numerically. They depend on Dg, 
Dp, I, and r. The range oi Dp we will consider is from 
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A to 3A. In Fig. 10, we plot these proximity lengths on 
both the F and S sides for three different values of /, at 
T = 4.0. The left panels show the /o proximity lengths 
and the right panels that extracted from /i . Recall that 
I — 1 corresponds to the half- metallic limit. We con- 
sider first the F side (top two panels) One can clearly see 
that both Lq^p and Li^p are correlated to the strength 
of the exchange fields. Fig. 10 displays a period of near 
A/2 for both Lo,f and Li,f at / = 1.0. We also see 
that the peak heights increase only slowly with increas- 
ing Dp. Also, the locations of the maxima or minima 
of Lo,F are locations of minima or maxima, respectively, 
of Li p. This is as one might expect from the rotating 
character of the field. On the other hand, for / = 0.1 
or / = 0.5 the periodicity is not clear since, for reasons 
already mentioned, the intermingling of periodicities be- 
comes more complicated. Overall the proximity lengths 
are larger than those in the half-metallic limit. However, 
one can still say that both Lo,f and Li^p gradually in- 
crease, although with fluctuations, with Dp. 

The superconductor on the S side is intrinsically s- 
wave but because of the F layer, triplet correlations can 
be induced in it, near the interface, as seen in Figs. 8 
and 9. Their decay, which is now monotonic, can be 



equally characterized by the proximity lengths defined in 
Eq. (11). Results are plotted in the bottom panels of 
Fig. 10. The minimum of io,s is, for all three values of 
I, at Dp = 23 which is near 2 A. The maxima of Lq.S 
for / = 0.1 arc at Dp = 16 and Dp = 28 which are 
not far from 1.5 A and 2. 5 A respectively. On the other 
hand, for Li s the maxima for / = 0.1 are at Dp — 19 
and Dp = 31 and there is a minimum at Dp = 26. The 
locations of these maxima are still near 1.5A and 2.5A 
and they are only slightly different than what they are 
for Lq^s case. If one recalls the above discussion of Fig. 3, 
maxima of Tc occur when Dp is close to an integer mul- 
tiple of A. Since a higher Tc is correlated with a higher 
singlet pair amplitude, this suggests again that there ex- 
ists a conversion between singlet and triplet Cooper pairs. 
The dependence of Lq.s and Li^s, at / = 1.0, on Dp is 
harder to characterize. This is because the high value of 
/ reduces the scale of the overall proximity effect in S (i.e. 
the depletion of the singlet amplitude). At / = 0.5, one 
still finds that the approximate periodicity of Lo,S and 
Li_s is about A/2. The proximity lengths Lq^s and Li^s 
are again anti-correlated at / = 0.5: the maxima (min- 
ima) locations of Lq.s are near the minimum (maximum) 
locations of -L1.5. 
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FIG. 11. (Color online) Triplet proximity lengths vs. Dp for 
Lo,F and -£/o,s, and the right panels show Li,f and Li.g. The 

Recent experiments^^ in systems that consist of two 
superconducting Nb electrodes coupled via a Ho/Co/Ho 
trilayer have revealed that the long range effect of triplet 
supercurrents was much more prominent at particular 
thicknesses of the Ho layers. The magnetic coherence 
length in Ho in the experiment was ^ 5 nm which would 
correspond in our notation to / ~ 0.1.^'^ In the experi- 
ment, the Ho thickness was symmetrically varied and the 
critical current, Ic, at T = 4.2X was measured. Peaks 
of Ic corresponding to Dp = 0.5A and Dp = 2.5A were 
found. These experimental findings are consistent with 
our theory. Here, we have shown (see Fig. 10) that Li,s 
has maxima near 1.5A and 2.5A when / = 0.1 in the Dp 
range we have considered. We found another maximum 
a.t Dp ^ 0.5 A, not included in the range shown. The 
penetration lengths associated with S are as important 
as those associated with F when discussing the triplet 
proximity effect, because the system can open up the 
corresponding channels only when both of them are long 
ranged. We believe that no obvious peak near 1.5 A was 
observed because of the layout of their symmetric sys- 
tem. Therefore, one can conclude that the spiral mag- 
netic structures play an important role in the triplet prox- 
imity effects. Both experiment and theory confirm that 
the existence of the long range proximity effects depends 
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I = 0.5, Ds = 150, at different t values. The left panels show 
es are guides to the eye. 

on the relation between the thickness of the magnetic 
layers and the wavelength of their magnetic structure. 

Having seen in the previous two figures that triplet am- 
plitudes may substantially pervade even rather thick Ho 
layers at moderate values of r, it is of interest to investi- 
gate the T dependence of the proximity lengths in these 
nanoscale F/ S systems for times roughly up to 2tt, in our 
dimensionless units. We therefore present in Fig. 11, the 
triplet proximity lengths as a function of Dp for / = 0.5, 
and at different values of r. The panel arrangement is as 
in the previous figure. Thus, in the top panels where we 
plot Lq p and Li p, we see that both of them depend only 
weakly on r, in the range considered. This is in part be- 
cause of the relatively thin F layers included in the plot. 
The triplet amplitudes vanish at t = but can saturate 
quickly through the F region as soon as r increases. In 
contrast, on the much thicker {Ds = 150) S side (bottom 
panels) both Lq.s and Li^s increase with r, as is consis- 
tent with expectation and previous work involving F/S 
systems with misaligned exchange fields. Furthermore, 
the overall shape of the proximity lengths on the S side 
does not change with r and only the magnitude evolves. 
Quite remarkably, the minima of Lq s and ii,s are very 
deep, and the value of these lengths at their minima is al- 
most r independent and nearly the same at all minima in 
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FIG. 12. (Color online) Normalized (see text) local magne- 
tization components plotted as a function of Y for several 
values of I. From top to bottom, x, y, and z-components are 
shown. We use Dp = A and Ds = 150 in this figure. 



the range plotted. The minima are separated by A/2. If 
one compares the left and right panels one can see that 
the locations of maxima in one approximately coincide 
with the position of minima in the other: the left and 
right panels are again complementary to each other as 
was the case with the plots in Fig. 10. 



D. Local magnetization and LDOS 

Next, we discuss other important physical quantities 
that are related to the proximity effects including the lo- 
cal magnetization, m{y), and the the local DOS (LDOS). 



Considering that the ferromagnetism can drastically alter 
the superconductivity, one might wonder about the oppo- 
site case: how the local magnetizations behave near the 
F/S interface. These so-called reverse proximity effects 
have been studied^'^^'^^'^^"^° for a number of multilayer 
F/S configurations with uniform exchange fields in each 
magnetic layer. Here, the space-varying exchange fields 
in F oscillate in the x — z plane and are constant along 
the y direction (see Fig. 1). We computed the local mag- 
netizations using Eq. (10) for Dg = 150, Dp = A and 
three different values of /. The results are normalized 
in the usual^^'^° way so that for a putative bulk F ma- 
terial with a uniform internal field characterized by the 
parameter / the quantity plotted would have the value 

[(1 + 1)3/2 _ (1 _ /)3/2]/[(l + J)3/2 j)3/2]^ j^^ 

Fig. 12, each component of m is shown in a separate 
panel and their behavior plotted throughout the whole 
spiral magnet region and some distance into S near the 
interface. Consider first the x component: The corre- 
sponding component of the internal field (see Eq. 1) van- 
ishes at the outer interface {Y = 0) and goes smoothly 
to zero at the F/S interface which in this case is at 
Y = Dp = A. As a consequence, one can see that the 
rrix component undergoes a full period of oscillation in 
the F material. The maximum and minimum values as a 
function of / arc numerically what they should be, given 
our normalization. However, as Fig. 12 clearly shows, 
the self-consistently determined does not vanish at 
the F/S interface and instead penetrates a short dis- 
tance inside 5*. This is a manifestation of the reverse 
proximity effect. For the other transverse component, 
rUz the situation is more complicated. The field compo- 
nent hz , out of phase with hx , does not vanish smoothly 
at y = nor at Y — Dp- Therefore the corresponding 
•niz component in F is squeezed, and in addition to the 
peak at Y = A/ 2, which has the expected location and 
value, there are two smaller peaks at intermediate val- 
ues. At the interface between materials, penetration of 
this component is appreciably more considerable than for 
rrix. The longitudinal component niy, which is induced 
by the uniform hy component, behaves qualitatively as 
transverse components do in uniform ferromagnet F/S 
structures. Penetration into the S layer occurs over a 
relatively short distance, except at the smallest value of 
/ where it is relatively larger although the overall scale is 
of course smaller. The value of my in the F layer is again 
the expected one, consistent with our normalization. 

Finally, we wish to discuss the LDOS. Here we will 
present results for the LDOS, as given in Eq. (9), summed 
over spins, integrated over cither the F or the S layer, 
and normalized, as usual, to its value in the normal state 
of bulk iS* material. The results are given in Fig. 13 where 
the energy scale of the horizontal axis is in units of the 
superconducting gap of bulk S material, Aq. The left 
panels of Fig. 13 show the LDOS integrated over the F 
(top) and S regions (bottom) for Dp = A, 1.5A, and 
2A. The superconductor has a thickness Dg ~ I.SSq 
and F has a relatively weak exchange field, / = 0.1. For 
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FIG. 13. (Color online) LDOS averaged over the F regions (top panels) and S regions (bottom panels), plotted vs. energy. 
On the left panels, the integrated LDOS is shown for different Dp and I = 0.1. On the right panels, the integrated LDOS is 
shown for different I and Dp — A. In all cases, the superconductor width is set to Ds = 150. 



Dp ~ A, one can clearly see, for the integrated DOS in 
the S side, peaks near e/Aq ~ ±1 as in the ordinary bulk 
spectrum. There is additional subgap structure including 
proximity induced bound states at smaller energies fol- 
lowed by a very deep dip-nearly a minigap. Overall, the 
DOS structure contains traces of the familiar DOS for a 
pure bulk superconductor. On the F side, the integrated 
LDOS at this value of Dp still exhibits BCS-like peaks 
at e/Ap = ±1 and subgap dip, but the whole structure is 
much weaker and the depth of the dip much smaller. It 
is indicative of the superconducting correlations present 
in the F region. In contrast, the subgap superconduct- 
ing features in the integrated LDOS for larger Dp values 
(Dp = 1.5A and 2A) are much less prominent, although 
the peaks near e/Aq = ±1 remain. Nonetheless, there 
are still shallow and discernible signatures in the gap re- 
gion, on both the F and 5* layers. For these two larger 
values oi Dp the results (as compared on the same side) 
are remarkably similar. This is surprising at first since 
we have already seen that Tc{Dp) has in this range of / 
maxima near Dp ^ A and Dp ^ 2 A and a local min- 
imum at Dp = 1.5A as shown in Fig. 3. From that, 
one might naively guess that the integrated LDOS for 
Dp = 2A should behave as that at Dp = A, with a dif- 



ferent integrated LDOS ior Dp ~ 1.5 A. This expectation 
is incorrect because, as one can see on a closer inspection 
of Fig. 3, Tc at D^F = A is higher than that at Dp ^2A 
although both are near local maxima. Furthermore, Tc 
a.t Dp ^ 1.5 A is closer to the value a.t Dp = 2A than 
to that at Dp = A. Since Tc is associated with the mag- 
nitude of the singlet pair amplitudes, in which the LDOS 
is indirectly correlated to, one should conclude that the 
LDOS corresponding to Dp = 2A should be similar to 
Dp = 1.5A rather than Dp = A. Indeed, the results 
confirm this notion. 

On the right panels of Fig. 13, we present the inte- 
grated LDOS on both the F and the S sides for differ- 
ent exchange fields, / = 0.1, / = 0.5, and / = 1.0, at 
Dp ~ 12 = A. We see that when / is increased from 
/ = 0.1, the integrated LDOS on the F side becomes 
quite fiat (at the value (1/2) [(1 + 7)^/2 + (1 - 7)^/2] as 
per our normalization) and essentially devoid of a su- 
perconducting signature. On the S side, the integrated 
LDOS at / = 0.5 and / = 1.0 still retains some vestiges 
of the structure seen in the / = 0.1 case. However, the 
integrated LDOS at / = 0.5 on the 5* side is slightly dif- 
ferent than that at / ~ 1.0. The dip in at / = 1.0 is 
wider than for / — 0.5, in a way more superconducting- 
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like. What happens is that at larger values of / the 
mismatch between the Fermi wavevector in S and the 
Fermi wavevectors in the up and down bands in F in- 
creases. This diminishes the penetration of the Cooper 
pairs into S and hence the overall scale of the proximity 
effects. We recall that the overall dimensionless scale of 
the proximity effect in F is roughly = 1//. Conse- 
quently, superconductivity is impaired in 5* over a smaller 
scale when it is in contact with a stronger ferromagnct. 
Having said that, one might argue that at / = 0.1 the 
integrated LDOS on the S side should have a smaller dip 
than the other two curves for stronger /. However, we 
have to consider here also the overall behavior of the Tc 
vs. / curves at constant Dp- This behavior is once again 
oscillatory but with a superimposed decay. The overall 
decay results in Tc being higher at / = 0.1 than at either 
/ = 0.5 or / = 1, but the oscillations produce a higher 
value of Tc at / = 1 than at / = 0.5. This explains the 
progression of the curves. All the above discussion and 
results indicate that the LDOS can provide, if properly 
analyzed, another perspective and additional information 
on the superconducting nature of our bilayers. 



IV. CONCLUSIONS 

We have studied several aspects of proximity effects in 
F/S bilayers, where the ferromagnet has a spiral struc- 
ture characteristic of rare earths such as Ho, by numeri- 
cally solving the self consistent BdG equations. We have 
calculated Tc{Dp), the critical temperature as a function 
of magnet thickness, for different parameter values. The 
Tc{Df) curves exhibit a fairly intricate oscillatory behav- 
ior which is found to be related to both the strength / 
(as they would for a uniform magnet) and the periodic- 
ity A of the spiral exchange fields inherent in the magnet. 
As is the case for F/S structures in which F is uniform, 
we observe reentrant behavior with Dp when / is strong 
enough. The physical reason behind this Dp recntrance 
in our bilayers is similar to that in ordinary F/S struc- 
tures but the additional periodicity associated with the 
magnet, which in many cases dominates the oscillations, 
makes the behavior more complicated. As a function of 
Ds, we find that Tc{Dp) can also exhibit Dp reentrance 
even at small / when Ds is of the order of the super- 
conducting coherence length. The additional oscillations 
produced by the magnetic structure lead also to pure 
reentrance with temperature: superconductivity occurs 
in a finite temperature range Td < T < Tc2- An example 
of this reentrance at a very small Dp (Dp ^ 0.5A) was 



previously^® presented. Here we report that this recn- 
trance can also occur when Dp > A, where it should be 
experimentally easier to realize. To elucidate the physics 
underlying these reentrant phenomena, we have eval- 
uated the singlet pair amplitudes and thermodynamic 
functions. The competition between condensation en- 
ergy and entropy is responsible for the T reentrance: 
the superconducting state may be, under certain circum- 
stances, the high entropy state, leading to recovery of the 
normal state as T is lowered. The calculated thermo- 
dynamic quantities are fully consistent with the Tc{Dp) 
phase diagrams and the singlet pair amplitudes. 

When the magnet has a spiral structure both the 
m = and m = ±1 odd triplet components can be in- 
duced simultaneously. This is not the case in uniform- 
magnet bilayers: at least two uniform misaligned F layers 
are needed to generate the m = ±1 component. We stud- 
ied the odd triplet pair amplitudes in our bilayers, and 
found them to be long-ranged in both the S and F lay- 
ers. We have analyzed the time delay dependence of the 
odd triplet amplitudes. The results are consistent with 
our previous work on both F1/S/F2 and F1/F2/S trilay- 
crs, but the additional A periodicity leads to important 
differences. We characterized the triplet long range be- 
havior by introducing the appropriately defined lengths. 
We found that the relevant proximity length oscillates 
with Dp and these oscillations depend on the strength 
and periodicity of the exchange field. Our methods are 
likely appropriate for many experimental conditions, as 
evidenced by the consistency of our results with recent 
tunneling experiments. '^'^ 

We have also considered the reverse proximity effects: 
the influence of the superconductivity on the magnetism. 
We found all three components of the local magnetization 
penetrate in slightly different ways into the S layer. At 
larger I this is a short-ranged phenomenon, but it is oth- 
erwise for weak magnetism. Both mx and oscillate in 
the F regions to reflect the spiral exchange field. Finally, 
the calculated LDOS reveals important information and 
discernible signatures linked to the proximity effects in 
these bilayers and are correlated to the superconducting 
transition temperatures. 
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